Impurity scattering and localization in d-wave superconductors 
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Strong evidence is presented for the localization of low energy quasiparticle states in disordered 
d-wave superconductors. Within the framework of the Bogoliubov-de Gennes (BdG) theory applied 
to the extended Hubbard model with a finite concentration of non-magnetic impurities, we carry out 
a fully self-consistent numerical diagonalization of the BdG equations on finite clusters containing 
up to 50 X 50 sites. Localized states are identified by probing their sensitivity to the boundary 
conditions and by analyzing the finite size dependence of inverse participation ratios. 
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In conventional s-wave superconductors nonmagnetic 
impurities have little effect on the superfluid density and 
the transition temperature, which may be understood 
from Anderson's theorem [Q. The situation is dramati- 
cally different in the high temperature cuprate supercon- 
ductors where nonmagnetic impurities exhibit a strong 
pair breaking effect 0. This is interpreted as a man- 
ifestation of unconventional, most likely d-wave pairing 
symmetry. Understanding the role that impurities play 
in such superconductors is crucial for the interpretation 
of experimental data. A good example of this is the inter- 
pretation of experiments measuring the temperature de- 
pendent penatration depth |^,^ in YBa2Cu307_:c- The 
behavior of various quantities in the presence of disorder 
also serves as an important test for theoretical models of 
microscopic pairing interaction. 

In the present paper we concentrate on the possibil- 
ity of localization of the low-energy quasiparticle states 
in d-wave superconductors with line nodes in the gap 
on the fermi surface. We wish to address the question of 
whether these low lying states are strongly localized with 
a relatively short localization length or whether they 
are essentially extended. This is a question of significant 
importance for many experiments since the transport 
and thermodynamic properties are largely determined by 
these low energy excitations. 

The problem of localization in d-wave superconductors 
was first considered by Lee Q who, appealing to argu- 
ments from the scaling theory of localization, found that 
in the limit of unitary (i.e. strong) scatterers the quasi- 
particle states are strongly localized below the mobility 
gap 7o, even if the impurity concentration is sufficiently 
small that the normal state wavefunctions are essentially 
extended. For moderate concentrations of strong scatter- 
ers, 7o is a reasonable fraction of the maximum gap which 
allows for the possibility of experimental confirmation. 
One consequence of such a scenario is the prediction of a 
universal limit of conductivity cr(w ^ 0) ~ (e'^ /2Trh)S,o/a, 
independent of the scattering rate (^o is a coherence 
length and a is the lattice constant). Infrared reflectance 
data on ion-irradiated YBa2Cu307_a; have been inter- 
preted using this picture 

More recently Balatsky and Salkola Q presented re- 
sults that contradict this scenario. Their argument is 



based on a physical picture of a single impurity wave- 
function which, according to an earlier calculation within 
the self-consistent T-matrix approximation is highly 
anisotropic with slowly decaying (~ 1/r) tails along 
the (11) and (11) diagonals. Overlaps between these 
tails lead to strong interactions between the quasiparticle 
states on distant impurities, which then form a network 
of extended impurity states capable of carrying current. 
This network percolates across the entire system and in- 
hibits the localization by disorder at the lowest energies, 
creating an "inverse mobility gap" 7c < 7o- 

Both theories described above resort to approxima- 
tions when treating the scattering from impurities and 
they ignore the gap relaxation near the impurity sites 
resulting from their pair breaking property. Since the 
arguments for and against localization are quite subtle 
it may well be that the above mentioned details are im- 
portant, in which case perhaps the best way to address 
the problem is through numerical calculations. Hatsugai 
and Lee Q studied localization numerically in a simpler 
but related model of Dirac fermions on the lattice, whose 
excitation spectrum is similar to that of a d-wave su- 
perconductor. By examining the sensitivity of the wave- 
function to boundary conditions they concluded that the 
low energy states are indeed strongly localized. This 
would seem to be in qualitative agreement with the the- 
ory of Lee , however since the system possessed no off- 
diagonal long range order, the quantitative predictions of 
this theory could not be verified. While this work was 
strongly suggestive, without treating a system with su- 
perconducting order one cannot establish the existence 
of localization in a d-wave superconductor. A convenient 
framework for studying the effect of impurities within 
a simple lattice model of d-wave superconductivity was 
outlined by Xiang and Wheatley The same model 
was later used to study localization induced by weak 
disorder, modeled by a white noise random component 
added to the chemical potential This corresponds 

to the scattering in the Born limit where the localiza- 
tion effect had been predicted to be neghgible The 
numerical work indeed confirmed that the localization 
length in the superconducting state remains comparable 
to its normal-state value (i.e., very large) down to the 
lowest energies. 
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In the present work we address the case of a di- 
lute density of strong (nearly unitary) scatterers. It is 
this limit that has precipitated controversy 1^,01, and 
that is most relevant experimentally (e.g., for Zn doped 
YBa2Cu306-5)- In order to confront the question of lo- 
calization in such systems we solve a simple tight bind- 
ing model with an on-site repulsion and nearest neighbor 
attraction which give rise to superconductivity in the d- 
wave channel, on finite clusters using the selfconsistent 
Bogoliubov-de Gennes (BdG) technique. The main ad- 
vantage of such an approach is that the impurity scat- 
tering is treated exactly and that it is directly relevant 
to c?-wave superconductors. Our principal result is that 
quasiparticle states are strongly localized at low energies 
below a mobility gap 70 . The value of 70 , as well as the 
localization length S,Li are in a good agreement with the 
theory of Lee Q . 

The Hamiltonian we consider has been used previ- 
ously to study the vortex structure [ p^|Jl^ and impurities 
| pO| , pl| in a d-wave superconductor: 

H ^ - tY^ c\^Cj„ - ^ n,^ -f ^ yr"^ni„ 

{ij)a icr ia 

-h Vb^n^tn^x -I- Y ^ niafija' (1) 

Here (ij) stands for nearest neighbor pairs, and the nota- 
tion is otherwise standard. The last term models strongly 
repulsive impurities: V^"^^ — y""P > at randomly cho- 
sen sites with density riimp and V^"^^ = on all other 
sites. If one takes Vb > and Vi < this model gives 
rise to pairing in the d-wave channel. In mean field the- 
ory (|^) can be solved by defining the pairing amplitudes 

Ao(r,) = Vo{c,ic^i), Asir,) = Vi(c,+5|C4), (2) 

where S = ii, ±y are nearest neighbor vectors for a 
square lattice. The resulting mean field Hamiltonian is 
then diagonalized using the Bogoliubov transformation 



n 
n 

to the quasiparticle operators jna- Within mean field 
theory the Hamitonian (|l]) is diagonalized when m„ and 
Vn satisfy the BdG equations p^ , p^ 

where 

iun{ri) = -i^u„(ri + S) + (V""'^ -^)u„(ri), 
s 

Au„(r,;) = Ao(r,) + ^ A5(r0w„(r, -I- S), (4) 

i5 




FIG. 1. Amplitudes of a) d-wave and b) s-wave order pa- 
rameters in the vicinity of a single impurity for a 17 x 17 
system. Note the different scales for d and s. 

subject to the constraints of self-consistency 
Ao(r) = VoY,Mr)vlir) tanh(S„/2fcBT), 

n 

A^(r) = Y E [""(^ + -^X^^) + Mr)v:{r + 6)] 

n 

xtanh(£:„/2fcBT), (5) 

where the summation is over positive eigenvalues En only. 

For a system of linear size L, solving the system of 
equations with periodic boundary conditions requires 
diagonalizing a 2L^ x 2L^ matrix. For a suitably chosen 
initial order parameter distribution we solve the system 
(^) using the standard LAPACK diagonalization routine. 
We than compute new gap functions from Eqs. (|^) and 
iterate until the desired convergence is achieved. Typ- 
ically 7-8 iterations are necessary to establish five-digit 
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convergence in the free energy and the average gap. Pe- 
riodically we perform longer runs of 25-30 iterations to 
confirm the accuracy of the solutions. No significant de- 
viations are found between short and longer runs. 

We have solved the BdG equations (|3|) at T = for 
the following set of model parameters: fi — —t, corre- 
sponding to the band filling factor n ~ 0.68 {n = 1 
is a half- filled band); Vq = 3t, and Vi = —A.bt. In 
the absence of disorder Eqs. (^) are easily solved ana- 
lytically by appealing to translational invariance. One 
obtains the usual BCS type excitation spectrum _Ek ~ 
(e^ -I- |Akp)^/^ with Ek — —2t{coskx + cosky) — jj, and 
Ak = 2Ad(cos kx cos/cy). For the above parameters we 
obtain A^ = 0.67t. With respect to real materials this 
value is exaggerated; however, our choice of parameters 
was motivated by the desire to have the large mobility 
gap 70 (which according to Lee Q scales with A^^), nec- 
essary to study localization numerically in a finite system. 

In Fig. we display the amplitudes of d- and extended 
s-wave order parameters defined as p^ ] 

d{r) = i[A,(r) + A_,(r) - Ay{r) - A_g(r)] 

sir) = ^[A£(r) + A_£(r) + Ay{r) + A^y{r)] 

in the vicinity of a single impurity with — lOOt. The 
d-wave order parameter is suppressed at the impurity site 
and recovers its bulk value over 2-3 lattice spacings, just 
as one would expect for a pair-breaking impurity. The 
s-wave, which vanishes in the bulk, is nucleated near the 
impurity site and displays an interesting real space struc- 
ture with a "d-wave" symmetry: it vanishes along the 
|a;| = |y| diagonals and changes sign upon 90° rotation. 
This behavior can be easily understood from the corre- 
sponding Ginzburg-Landau theory We note that 
being of the form d ± s this state does not break time 
reversal symmetry. 

For a finite density of impurities we probe for local- 
ization in two different ways. First, for a given system 
size and disorder configuration we evalute a generalized 
inverse participation ratio |I6{, 



((KP) + (KP))^ 



(6) 



where (...) stands for the sum over all sites. As a func- 
tion of increasing system size L this quantity decreases 
as ~ 1 /L^ for extended states and approaches a constant 
value ~ (a/Ci)^ for ^ state localized within the char- 
acteristic length Second, we study the sensitivity 
of states to boundary conditions by applying a uniform 
phase twist x accross the length of the system. We con- 
sider wavefunctions with the "twisted" boundary condi- 
tion along the a;-direction, ^'(a; -|- L, ?/) = e'^^(a:,j/), and 
periodic along y, and we compute the stiffness |p"6|j9| of 
the n-th eigenvalue 
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FIG. 2. Inverse participation ratio a„ and stiffness k„ 
plotted as a function of energy for L = 40 and Ui = 0.06. 
Dashed line marks the mobility gap 70 = 0.84 estimated from 
the theory of Lee j^. Inset: Finite size scaling of the inverse 
participation ratio averaged over the energies below (a^) and 
above (a^) the mobility gap 70. The error bars reflect the 
scatter of the data from 6 different impurity conflgurations. 
Solid and dashed lines are guides to the eye only. 



Extended states will be sensitive to twist, and thus k„ is 
expected to be large. Localized states with <C i will 
be insensitive to twist and k„ is expected to be small. 

We have studied systems with = 0.015 and rii — 
0.06 of strong repulsive impurities with = 100<. In 

Fig. 1^ we display typical results for the inverse partici- 
pation ratio a„ and stiffness function of energy 
for rii = 0.06. There is a pronounced qualitative differ- 
ence in a„ between the low energy states below E/t c^i 1 
and the high energy states. This is strongly suggestive 
of a mobility gap with a magnitude close to that pre- 
dicted by Lee; [|| 70 ~ 0.84t for the parameters listed 
above. While not so pronounced, this mobility gap can 
also be observed in the behavior of k„. For lower density 
of impurities, Ui — 0.015, we obtain similar results with 
the apparent mobility gap reduced by roughly a factor 

of two. This is in agreement with Lee's prediction that 

1/2 

70 - n/ . 

Even stronger evidence for the existence of the mo- 
bility gap can be obtained by analyzing the system size 
dependence of a„. To this end we have carried out cal- 
culations for L = 20,24,28,32,36,40,50 and six inde- 
pendent impurity configurations for each size. The inset 
of Fig. ^ shows two quantities, a"^ = {{an)) ek-^o 

= {{ayi)) E>-ia^ defined as averages of the inverse par- 
ticipation ratios a„ taken over the indicated range of en- 
ergy, as a function of I/L^. Data points for lie on 
a straight line which extrapolates to zero for L 00, 
just as one would expect for the extended states. Val- 
ues of appear to extrapolate to a finite value of 
~ 0.003, which implies that these states are local- 
ized with ^l/o ~ 20. We note that this value of is in 
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a reasonable agreement with the rough estimate given in 
Ref. [4], which gives ^l/o « 38. 

We have inspected visually the amplitudes of wave- 
functions u,i(r) and v„(r) for signs of localization. We 
have found that, in agreement with the above analysis, 
these states are spatially localized below the mobility gap 
with the characteristic length scale of about 20 — 30a, and 
appear to be extended at higher energies. 

The question naturally arises why we do not observe 
delocalization effects below the inverse mobility gap 7c 
predicted in Ref. [Q . In order to clarify this issue we have 
studied quasiparticle states associated with a single im- 
purity. For the model parameters described above {i.e., 
strong coupling) we found no evidence for the anisotropic 
cross-shaped states predicted in Ref. Thus, it is not 
surprising that for a finite density of impurities the de- 
localization mechanism based on long-ranged overlaps is 
inoperative and all we find is a sharp mobility edge below 
which all the states are strongly localized. Furthermore 
we found that in this regime it is very important to treat 
the order parameter self-consistently; a separate calcula- 
tion with spatialy uniform (i.e., not relaxed) order pa- 
rameter yielded wavefunctions with different spatial dis- 
tributions and energy eigenvalues. 

For the sake of completeness the above analysis was 
repeated for a different set of model parameters which 
may be closer to the realistic values for cuprates p3| : 
Vo = l.OSt, Vi = -1.05t, fi = -O.m, which imply 
smaller magnitude of the bulk gap. Ad ~ 0.066t. In 
this case we have indeed found the predicteded highly 
anisotropic quasiparticle state, however only the electron 
part of the wavefunction, Vnir), exhbited the predicted 
tails along the diagonal directions. The amplitude of the 
hole part, zin(r), was found to vanish along these diago- 
nals with two tails running parallel to it. Such a structure 
was recently also obtained from the T-matrix approxima- 
tion and it is not clear how it affects the analysis of 
the localization problem in a system with finite impurity 
concentration. Within our numerical method we were 
unable to address this question since for the above pa- 
rameters the estimated localization length j|] would be 
much larger than the maximum system size we can han- 
dle. Since in this case virtually identical results were 
obtained with a spatially uniform ordcrparameter, the 
method employed previously by Xiang |11| may be suit- 
able to address this problem. 

Thus, the evidence presently available allows for two 
distinct scenarios with regards to the existence of the in- 
verse mobility edge. Either (i) there is a crossover from 
the strong coupling regime which has (as we have explic- 
itly demonstrated) 7c = to the weak coupling regime 
with 7c > 0, or, (ii), it may be that 7c = for all coupling 
strengths which could be explained as a consequence of 
the complicated structure of the quasiparticle states dis- 
cussed above that was not included in the original anal- 
ysis of the many impurity problem 0| . 

To conclude, we have presented clear evidence for 
strong localization by impurities in d-wave supercon- 



ductors at low energies. Compared to earlier analyti- 
cal work 1^,0 our approach treats impurity scattering 
exactly within BdG theory and accounts for order pa- 
rameter relaxation near the impurity sites in a fully self- 
consistent way. The present model can easily be extended 
to finite temperatures and frequencies, thus opening posi- 
bilities for the analysis of the temperature dependent 
penetration depth, infrared conductivity, critical temper- 
ature, and other physically interesting quantities in the 
presence of disorder. 
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